Hysteretic dynamics of domain walls at finite temperatures 
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Theory of domain wall motion in a random medium is extended to the case when the driving field 
is below the zero-temperature depinning threshold and the creep of the domain wall is induced by 
thermal fluctuations. Subject to an ac drive, the domain wall starts to move when the driving force 
exceeds an effective threshold which is temperature and frequency-dependent. Similarly to the case 
of zero-temperature, the hysteresis loop displays three dynamical phase transitions at increasing ac 
field amplitude Hq. The phase diagram in the 3-d space of temperature, driving force amplitude and 
frequency is investigated. 
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Pinning dominated driven dynamics of elastic media 
in random environment is a paradigm for a vast diversity 
of physical systems. Examples include vortices in type II 
superconductors, charge density waves (CDW) in solids, 
stripe phases, Wigner crystals, dislocations in crystals, 
domain walls in magnets and many others [Q. Having 
appeared first in the context of dislocation dynamics , 
the scaling theory of glassy dynamic state of random elas- 
tic media came to fruition in the context of CDW ||] and 
vortex lattices in high temperature superconductors [^,^ , 
and enjoyed an impressive success in explaining a wealth 
of phenomenology of the low temperature vortex state 
1^. A closely related subject is the zero temperature 
depinning transition first studied for CDWs and do- 
main walls [f|j9|. Despite the significant recent progress, 
several key questions specific to glassy dynamics are yet 
poorly understood. One of such fundamental key issues, 
although known and extensively studied for more than 
hundred years in magnets, is hysteresis of interfaces sub- 
ject to the applied ac drive and related aging and memory 
effects. A quest for urgent progress in understanding hys- 
teretic behavior of magnetic domain walls is motivated 
also by emerging technological nano-scale magnetic sys- 
tems whose ac properties are controlled by the hysteretic 
dynamics of interfaces. 

A step towards theoretical description of hysteretic be- 
havior of disordered interfaces has been undertaken in 
pof , where the cyclic motion of the domain wall at zero 
temperature under the ac field was investigated and the 
resulting magnetization hysteretic loop was described. A 
finite temperature may change drastically the interface 
dynamics: thermally activated creep motion becomes 
possible at any small drive. 

In this Letter we develop a unified description of ther- 
mally activated and over-threshold domain wall dynam- 
ics in impure magnets. We demonstrate that at finite 
temperature new scales of length, activation energy and 
force appear leading to emergence of a new, temperature 



and frequency-dependent threshold field in the case of 
ac drive. The latter is the first in a series of dynamical 
phase transitions. To be specific, we will speak on mag- 
netic domain walls. Accordingly we will be using either 
of terms "force" or "field" equivalently. 

Finite temperature dc dynamics The essential of the 
zero-temperature dynamic behavior of an elastic medium 
in a random environment is the existence of the finite 
threshold depinning force hp, separating immobile at 
h < hp and sliding at h > hp states of the system. Near 
the threshold the sliding velocity v shows a critical be- 
havior |^|-^ V ^ {h — hpY . At finite temperatures and 
h <^ hp thermally activated drift motion controlled by 
the static rugged energy landscape occurs. The latter is 
governed by the interface free energy 



n = J d^x!^^r{\/zf + v{x, z)-hz{x)^ 



(1) 



where F is the interface stiffness, h is the external driving 
force and ^(x, Z) is the random impurity potential. D- 
dimensional vector x is the coordinate along the inter- 
face, and Z is the coordinate of the transverse interface 
displacement. In the following we assume that the dis- 
order average of the random potential vanishes. There 
are two different types of impurities, random bond (RB) 
and random field (RF) type in terms of magnetic models. 
The RB potential obeys the Gaussian statistics with: 



Vrb(x, Z) VRsiO, 0) = <5(x) S{Z) 



(2) 



where = VqC. vq, c and I denote the strength, the con- 
centration and the correlation length of the impurity po- 

z 

tential. In the RF case Vrf{x, Z) = J h{x, Z') dZ' where 



the RF /i(x, Z) has properties similar to yRB(x, Z). 

The static interface in a random environment becomes 
rough. Its roughness obeys the scaling law p^ ]: 
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for RF and 



where the roughness exponent is C = ^-3^ 
C ~ 0.2083(4- D) for 4 - D < 1 and C = 2/3 for D^2, 
for RB impurities, respectively |14 13]. The rough con- 



figuration devefops over length scales L > Lp, where the 
Larkin length Lp is a distance at which a typical fluc- 
tuation of the pinning forces, balanced by elastic forces, 
produces the transverse displacement w ^ I: 



Lr, 



nr/v) 



2/(4-D) 



(4) 



Note that the typical slopes w{L) /L of the wall vanish for 
L ^ Lp since C < 1- The energy barriers which must be 
overcome to depin a segment of the wall with the linear 
size L is [01: 



EbAL)~Tp{L/LpY- x = D-2 + 2C. 



(5) 



Here Tp = EB,o{Lp) 



TPL^~^ 



is a typical pinning en- 



ergy on a scale Lp. At temperature T > Tp the effective 
force necessary for depinning drops rapidly with the tem- 
perature j2| . If an external driving force h is applied, the 
total energy barrier Eb{L) becomes 



EBiL)^EBAL)-hL°w{L) 

« EbAL) (1 - {L/L^f~i) 



(6) 



where Lh ~ Lp{hp/hy^^^ The energy barrier reaches 
its maximum EB,max at L w {j^^y^'^^^'^^'Lh which gives 

EB,rnax ~ Tp{hp/ hY and ^ = x/(2 - 0- 

For completeness we note here that equilibrium RF 
systems inc? = Z3-|-l = 2 bulk dimensions have no long 
range order at length scales L ^ w Zexpc(r/u)^/'^ 
pGf . For weak disorder S^np is very large and can easily 
exceed the system size. As soon as an external field is 
applied such that Lh < domain wall motion is dom- 
inated by the forces on scales L^ and the absence of true 
long range order can be ignored. 

As first found by Middleton jl^ , close to hp the small- 
est energy barrier vanishes as {h — hpY with 9 — 2. We 
therefore get the effective energy barrier at h hp if we 
replace in the above expressions Tp by fp = Tp{ ^^^ ^ )^. 

The time scale to overcome this barrier is of the order 
t{Lii) « To exp {j^B^max/T^, where Tq is a microscopic 

hopping time, which leads to an average velocity: 



v{h) « 7/1 exp 



T 



hp — h 



h 



(7) 



where 7 is the effective friction coefficient PP,p^. 

At low temperature T <^ Tp the dynamic threshold hp 
separates the creep regime from the active sliding regime. 
As can be seen from ^ a characteristic crossover field 
hx plays the role of the depinning force, where 



At h K hx, the drift velocity increases rapidly and at 
larger fields it displays almost linear behavior v « 7/1. 
Note that hx is a monotonously decreasing function of 
temperature with a maximum hx — hp at T = 0. 

In a close vicinity of the threshold field hp, the effec- 
tive energy barrier becomes small and even small thermal 
fluctuations may be sufficient to overcome it. At finite 
temperatures and h ^ hp the wall moves via thermal ac- 
tivation process with velocity given by eq. (|^). Strictly 
speaking, it means that at finite temperatures the critical 
point shifts from h — hp to h — 0. Yet there remains a 
memory of the critical behavior around h k, hp displaying 
itself in a crossover behavior at finite but low tempera- 
tures. The crossover is seen as a rounding of the h — v 
characteristics v{h hp) ~ T^^^. We now can write an 
interpolation formula for the velocity which is valid in a 
wide range of variables: 



v{h,T) = jhF{x,y); x = hjhp, y = Tp/T, 



F(x,y) = 



9(1 -x) 
1 + {yx-i'Y/<> 

+Q{x-l) 



exp 



-yx ^ {1 — x) 



(9) 



(10) 



l + [yx-t'Y/<> 
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Here Q{x) is the step function equal to zero at x < 
and equal to 1 at positive x. The interpolation for- 
mula ( p^ ) satisfies following requirements: (i) v{h, T) = 
7/1 at any fixed T and h ^ hx; (ii) v{h,T) = 



7/1 exp 



~ T \ hp 



for hp — h hp and T <C Tp] 



(iii) v{h,T) K, 7/1 exp 



T 



for T '^Tp,h'^hp 
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(v) 



and EB,max/T > 1 ; (iv) v{h,T) « jhp{h/hp - 
for {h/hp - 1) < 1 and T < Tp{h/hp - if; 
v{hp,T) « -fhp{T/Tp)l3/e for T « Tp. 

So far we assumed that the propagating interface is 
self-affine. This is confirmed by numerical simulations in 
D > 1 interface dimensions for systems with weak disor- 
der ||l8| . In = 1 dimensions the situation is less trans- 
parent: in simulations which use a bounded distribution 
of random fields the interface appears to be self-affine 
Jigf or faceted j2^ depending on whether lattice effects 
are avoided or admitted, respectively. We ignore here the 
possibility of faceted growth which occurs only in systems 
with narrow magnetic domain walls. For an unbounded 
distributions of random fields however a percolative self- 
similar domain wall propagation was observed pl[ |. In 
the following we will always assume, that the random 
fields distribution is bounded such that the domain walls 



remain well defined. This is also confirmed by our earlier 
simulation outside of the critical region . 

Alternating fields. If the external drive is oscillating 
with frequency lu, h = /igsinwi, the barriers for which 
Ltj t{L) > 1 cannot be overcome during one cycle of the ac 
field. From the condition cj r = 1 we find a new frequency 
and temperature dependent magnetic field hi^ which obeys 



hp \ T A hp 



(11) 



where A — lnl/(cjro). h^ plays the role of the dynamic 
threshold. At low fields ho < h^ there is no macroscopic 
motion of the wall, its segments oscillate between the 
metastable states with close energies giving rise to dis- 
sipation Drift of the wall starts at h^ > h^i- We 
assume lutq ^ 1, so that h^ < hr- Various regimes of 
domain wall motion are summarized in Fig.l. 
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FIG. 1. The phase diagram for the domain wall hystere- 
sis at finite temperature and frequency. Solid and vertical 
lines separate no sliding, thermal creep and mechanical drift 
regimes. See also explanations in the text. 

Having derived the domain wall velocity as a function 
of the driving field in the different h — T-regions we con- 
sider now the magnetic hysteresis following from the mo- 
tion of a single wall under the infiuence of an oscillating 
field. Since substantial length scales are larger than Lp, 
where slopes are small, the domain wall will be considered 
as a straight line (plane) characterized by one coordinate 
Z |10|] . Its dynamics is determined by equation of motion 



Z = v{h{t)) 



(12) 



Z varies between limiting values and L. Here L is the 
linear size of the sample in the case of a single domain 
wall or, in the multi-domain case, equal to the average 
distance between expanding nuclei . For harmonically os- 
cillating field h = ho sin uit, equation ( |l^ can be rewrit- 
ten in terms of h only: 



dZ 
dh 



vih) 



UJyJJl^~^h?' 



Z{h = 0) = 



(13) 



Equations ( p^ ) and ( p^ ) are valid for h > h^. The field 
region h < h^ where the motion has zero drift velocity 



will not be considered here. The value /ic^ plays the same 
role as the threshold field hp plays at zero temperature. 
At h > h,^ the hysteresis is dominated by the activation 
processes. In this respect it is similar to the nucleation 
dominated hysteresis described in ||l^ with two essential 
differences. First, the activation relates to the formation 
of a nucleus on the interface, not in the bulk. Second, 
the activation energy depends on the magnetic field as a 
power function /i^^ due to the distribution of barriers 
depending on their length scale. 
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FIG. 2. Schematic pictures of hysteresis loops(HL). (a) In- 
complete HL for ho < hti- (b) Symmetric HL for ho — hti- 
(c) The HL for hti < ho < ht2. (d) The HL for ho > ht2. 

h,. and ho are all marked in all figures 



The values hp , h 



We start with large temperatures T > Tp at which the 
zero-temperature threshold field hp plays no role (field 
sweep (2) in Fig.l). In this case the motion of the do- 
main wall is determined by equation ( 1^) for h > h^^. 
As in the case of zero temperature |10 three dynamic 
phase transitions take place when the amplitude ho in- 
creases gradually at a fixed value of frequency. These 
transitions change the shape (symmetry) of the hystere- 
sis loop. At the first of them proceeding at ho = h^ 
the hysteresis loop first appears, at smaller amplitudes 
ho < the magnetization remains unchanged. The 
hysteresis loop appearing at ho > h^^ is characterized by 
incomplete magnetization reversal and refiection symme- 
try h — > —h, M — > M as shown in Fig. 2a. This symme- 
try as well as incomplete magnetization reversal persists 
until the next dynamic phase transition at ft-o = ft-ti. 

At /i > hti the magnetization reversal becomes com- 
plete and hysteresis loop symmetry changes to inversion 
h ^ -h, M ^ -M (see Figs. 2b, 2c). The value hti is 
determined by a requirement that the domain wall pro- 
ceeds from one sample boundary to another for half a 
period. At the next dynamic phase transition the sym- 
metry of the hysteresis loop remains unchanged, but the 
part of the cycle becomes reversible. Visually the hys- 
teresis loop acquires characteristic "whiskers" as shown 
in Fig. 2d. The point of this transition ht2 is determined 
by a requirement that the domain wall proceeds from one 
sample boundary to another for quarter period. 
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Starting from the transition amplitude = hti each 
hysteresis loop goes through three important points. One 
of them is hun at which the motion of domain wall starts. 
Two others the so-called coercive field he and reversal 
field hr- At coercive field the magnetization turns into 
zero, at reversal field the magnetization becomes com- 
pletely reversed (see Fig. 2b-2d). Note that for /iq = hti 
hr = and he = ho; for ho = ht2 hr = Hq. All these 
fields can be found in our case. Equations for hti, ^t2 
are: 



htn 



v{h)dh 



ULUL 



1,2 



(14) 



For the case T > Tp the equations ( |l4| ) read: 

X 

riLjL 



2^hT 



g(x)^ / ye-y~\x^ -y^r^'^dx (15) 



h K, hp retains features of the critical behavior and ex- 
hibits the rounding of the h — v characteristic, according 
to ||l^. At finite frequencies uj, a new characteristic field 
hu < hp comes into play, and the transition from the 
sliding regime to pinning dominated activation motion is 
shifted to h^^. 
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where a;„ = htn/hr, n = 1,2; x^ = h^/hr = l/JS}/^. 
Thus, xi, X2 are functions of a dimensionless parameter 
u — ujL/^hT, where L is the size of the system or an av- 
erage size of domains. Its asymptotic at small u results in 
Xn ~ [ln(2/nu)]^^/^. The fields hti; ^t2 are close in this 
case: {ht2 ~ hti)/hti ~ ln2/(/xlnM). The requirement 
hti > is satisfied if ''yhxTo < L. The coercive field 
he and the reversal field hr are determined by equations: 
M(/i,) = 0; I M{hr) \=M,. 

A simple analysis at small u results in he ~ hti; hr ~ 
ht2- The area A of the hysteresis loop at u << 1 and 
ho > hti does not depend strongly on the amplitude ho 
and becomes A « AhrMs- Fig. 2 shows typical hysteresis 
loops and illustrates the geometrical meaning of the field 
hi^,hc,hr. The dependence of magnetization on magnetic 



- 1 



field is given by equation: M{h) = Ms ^Hihl 

Finally, in the range of moderately low temperature 
T < Tp the more complete expression (10) and (11) have 
to be used in integrating equation (14). 

It interesting to note that the similar (dynamic) tran- 
sition from incomplete to complete hysteresis was ob- 
served even in a standard simplistic mean-field model for 
pure magnets with the reaction described by the Bril- 
loin function [|2|. This suggests that this kind of dy- 
namic transition we discuss, may be a generic property 
of nonlinear systems. Another note is in order: hysteresis 
in large multidomain magnet samples is a very complex 
phenomenon and cannot be always reduced to motion of 
a single domain wall (see for example numerical simula- 
tions of random Ising model in [Til ] , where hysteretic and 
memory elects unlikely reducible to motion of a sible DW 
were revealed). 

In conclusion, we have investigated critical creep mo- 
tion at low, T ^ Tp, and at high, T > Tp, temperatures, 
Tp being the depinning temperature, and constructed the 
dynamic phase diagram. At low temperatures creep at 



[1] for recent reviews see e.g. M. Kardar, Phys. Rep. 301, 

85 (1998), D.S. Fisher , Phys. Rep. 301, 113 (1998), 
[2] L. B. loffe and V.M. Vinokur, J. Phys. C 20, 6149 (1987). 
[3] T. Nattermann, Phys. Rev. Lett. 64, 2454 (1990). 
[4] M.V. Feigelman, D.B. Geshkenbein, A.I. Larkin, and 

V.M. Vinokur, Phys. Rev. Lett. 63, 2303 (1989). 
[5] G. Blatter, M.V. Feigelman, V.B. Geshkenbein, A.I. 

Larkin, and V.M. Vinokur, Rev. Mod. Phys. 66, 1125 

(1998); 

[6] D.S. Fisher, Phys. Rev. Lett. 50, 1486 (1983). 

[7] O. Narayan and D. S. Fisher, Phys. Rev. B 46, 11520 

(1992) . 

[8] T. Nattermann, S. Stepanow, L.H. Tang, H. Leschhorn, 

J. Physique 2,1483 (1992). 
[9] O. Narayan and D.S. Fisher, Phys. Rev. B 48, 7030 

(1993) . 

[10] I.F. Lyuksyutov, T. Nattermann, and V.L. Pokrovsky, 

Phys. Rev. B 59, 2262 (1999). 
[11] J. P. Sethna, K. Dahmen, S. Kartha, J. A. Krumhansl, 

B.W. Roberts and J. D. Shore, Phys. Rev. Lett. 70, 3347 

(1993); K. Dahmen and J. P. Sethna, Phys. Rev. Lett. 

71, 3222 (1993) . 
[12] B. K. Chakrabarti and M. Acharyya, Rev. Mod. Phys. 

71, 847 (1999). 

[13] T. Halpin-Healy and Y-C. Zhang, Phys. Rep. 254, 215 
(1995). 

[14] D. S. Fisher, Phys. Rev. Lett. 56, 1964 (1986). 

[15] For a RG-derivation of |^ see L. Radzihovsky, APS March 

Meeting 1998 and P. Chauve, T. Giamarchi and P. Le 

Doussal, Europhys. Lett. 44, 110 (1998). 
[16] J. Villain, J. Phys. Lett. (Paris) 43, L551 (1982); G. 

Grinstein and S.-K. Ma, Phys. Rev. Lett. 49, 685 (1982). 
[17] A. Middleton, Phys. Rev. Lett. 68, 670 (1992). 
[18] L.Roters, A. Hucht, S. Lubeck, U. Nowak, and K.D. 

Usadel, Phys. Rev. E 60, 5202 (1999), B. Koiller and 

M.O. Robbins, arXiv:cond-mat / 000418; and references 

therein. 

[19] U. Nowak and K. D. Usadel, Europhys. Lett. 44, 634 
(1998). 



4 



[20] H. Ji and M.O. Robbins, Phys. Rev. B44, 2538 (1991). B. Drossel and K. Dahmen, Eur. Phys. J. B3, 485 (1998). 

[21] H. Ji and M.O. Robbins, Phys. Rev. B46, 14519 (1992), 



5 



